Dynamics of the Peierls-active phonon modes in CuGeC>3 
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We reconsider the Cross and Fischer approach to spin- 
Peierls transitions. We show that a soft phonon occurs only if 
Qo < 2.2 Tsp. For CuGe03 this condition is not fulfilled and 
the calculated temperature dependence of the Peierls-active 
phonon modes is in excellent agreement with experiment. A 
central peak of a width ~ 0.2 meV is predicted at Tsp. Good 
agreement is found between theory and experiment for the 
pretransitional Peierls-fluctuations. Finally, we consider the 
problem of quantum criticality in CuGeOa . 

PACS numbers: 67.57.Lm, 75.10.Jm, 75.25,+z, 
63.70.+h, 64.60.Cn, 68.35.Rh 

Introduction Structural phase transitions come essen- 
tially in two varieties, those with a soft phonon mode 
and those without phonon softening and a central peak 
0. Typically one associates them to displacive and to 
order-disorder transitions respectively, even though there 
is no strict formal distinction between displacive and 
order-disorder transitions. It came then as a surprise 
that the spin-Peierls transition in CuGe03 [||, which 
had been shown to be displacive j^], shows no phonon 
softening [|H]- Even worse, the Peierls-active phonon 
modes harden by about 5% — 6% with decreasing temper- 
ature U . It has been generally assumed, up to now, that 
this behaviour is inconsistent with the RPA-approach by 
Cross and Fischer (CF) to the spin-Peierls transition Q. 
Here we show, that the CF-theory is actually fully con- 
sistent with the experimental results for CuGeC>3. We 
show, that soft phonons occur within RPA only if the 
bare phonon frequency Oo satisfies fio < 2.2Tsp. For 
larger phonon frequencies the phonon does not soften 
and a central peak develops at the spin-Peierls transi- 
tion temperature, T$p- We then test the applicability 
of RPA to CuGe03 by calculating the pretransitional 
Peierls-fluctuations. We find good agreement with ex- 
periment. Finally, we note that a key ingredient of the 
CF-approach, quantum criticality, can be tested for in 
CuGe0 3 . 

RPA approach The retarded phonon Green's function 
D q (u) is given by 



D q (u) 



uj 2 - n$(q) - 2n (q)P q (L>y 



(1) 



were fl (q) is the frequency of the bare phonon with mo- 
mentum q. In RPA one approximates the phonon self- 
energy Pq(u>) by g^Xq^), where X<?( w ) i s the dynamical 
energy-energy correlation function and where g q is the 
electron-phonon coupling constant, given by 
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where we have used 
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(3) 



for the spin-phonon coupling within a linear-chain model. 
Here S„ are the spin operators at site n, u n the displace- 
ment operators for the normal coordinates of the Peierls- 
active phonon-mode |^,|| and M is the effective mass of 
the normal mode. 
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FIG. 1. Plots of Rex^(w) (solid lines) for T = T SP and 
T — 2Tsp, as a function of lu/Tsp- The temperature in- 
dependent LHS of Eq. (|), (to 2 - Q.l)/(2Q gl), is plotted for 
Qo/Tsp = 1.3, 2.6, 4 (dashed lines), with g n given by Eq. (Q). 
The filled circles denote the position of the phonon frequen- 
cies, u) n . The arrows indicate the shift of ui^ with decreasing 
temperature. Inset: ImxTr(w). 

At the spin-Peierls transition a spontaneous dimeriza- 
tion occurs below T$p, at q = n/c. In the following we 
set the lattice constant c to unity in the theory formulas. 
Cross and Fischer observed, that the correct functional 
from for x g (o>) (in the limit to — > 0) can be obtained from 
bosonization ||, 



Txq(to) = -2dh 
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where d « 0.37 is a constant depending weakly on the 
momentum cut-off, A = v s \q — tt\ is the lower edge of the 



1 



two-spinon continuum (v s is the renormalized spin-wave 
velocity), and 



h{k) 



1 

2^ 



e^sinhOr))- 1 / 2 . 



Tx,j=7r(w) is scale-invariant and a function of lu/(2ttT) 
only (independent of the spin-spin coupling J) . This be- 
haviour is characteristic of quantum critical systems ||. 
For any temperature T > we can expand Txir{oS) in 
lu/(2ttT) as 

T X -(_0 = - Xo - i X x (^) + X2 2 + • • ■ , (5) 

with xo ~ 0.26, xi ~ 0-81 and %2 ~ 2.2. The position of 
the poles of D n {Lo) are then determined by the roots 
of 

= J Re X7rM «-_ + _( k _j , (6) 



2^0.9 



where Slo = ^o( 7r )- Typical plots of the left and right 
hand side of (g) are presented in Fig. |l| 
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FIG. 2. The temperature dependence of the phonon fre- 
quencies lutt for various values of tto/Tsp. Qq is the bare 
phonon frequency. The arrows indicate the respective mini- 
mal phonon frequency. 

A spontaneous lattice dimerisation, i.e. a macroscopic 
occupation of the Peierls-active phonon mode, occurs at 
Tsp when (0) has a solution for u> = 0. This determines 
the transition temperature as Jq] 



2^ 

T'sp = 7^— Xo- 



(7) 



Remarkably, Eq. (0) is independent of J, due to the scale- 
invariance of T Xw (w). We compare the prefactor of the 
terms ~ uj 2 of the rhs and lhs of Eq. (0) and find that for 



l/(2 5 2 tt ) > X2 /(4ir 2 Tl P ), 



(8) 



Eq. (||) has a single solution for T — T$p and inspection 
of the temperature dependence of this solution for T > 
Tsp (compare Fig. (|])) shows that this root continuously 
connects to the T = oo solution, Iuilt—hx — — ^o( 7r )- 
In the parameter regime defined by Eq. (||) the phonon 
softens completely. We can use Eq. ([_) to eliminate 
from Eq. (ph. We obtain 



T S p>^ x f^-~0A6n , 
2?r V Xo 



f!o< 2.2 T SP , (9) 



for the soft-phonon regime. For Slo > 2.2 Tsp the Peierls- 
active phonon does not soften completely and may even 
become harder with decreasing temperature, as illus- 
trated in Fig. |[ Near T = T$p an additional central peak 
shows up, leading to the phase transition. For CuGeC>3 
there are two Peierls-active phonon modes with energies 
|| _>i = 151 K and u>2 = 317 K respectively. Since (see 
below) f2 7 w w 7 (7 = 1,2) and Tsp = 14.1 K we find 
that CuGeC>3 is in the central-peak regime. 

Application to CuGeO^ In order to compare more in 
detail with the experimental results for CuGeOa we have 
generalized Eq. ([!]) for the case of two phonon frequen- 
cies. Denoting by D 1 / 2 (oj) the retarded Green's functions 
of the first and second phonon with bare frequencies fii 
and SI2 respectively, and by g\ and gi the respective spin- 
phonon coupling constants, we obtain 



D 1 (w)=I>i 0) (_.) + 



(.9 2 M 0) H+.9 2 2 ^ 0) H)x^M 



and an equivalent equation for D-z(ui). Here D^ 2 (lu) = 
2Q,i/2/ '{u 2 ~ ^1/2) @- An analysis similar to the one- 
phonon case can be performed for D\(uj) + Z)a(w). One 
finds 



Tsp= U~ + t^'* 

for the transition temperature and 

Ix^O^ g^+g^Ox 



T S p > 



2tt Vg 2 3 2 +g 2 2 0l 



(10) 



(11) 



for the soft-phonon regime. 

In order to determine g\ and 52 for CuGeOa we note 
that the lower/upper phonon mode contribute to the 
structural distortion below Tsp with weighting factors 
2 and 3 respectively || . This leads to 
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Eq. ( |l2] ) and Eq. ( ]To| ) determine the spin-phonon cou- 
plings 5 i, g 2 . For T S p = 14.1 K we find Oi = 3.15 THz 
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and Q2 = 6.61 THz for the bare phonon frequencies and 
91 = 0.86 THz (32 = V^gi) for the spin-phonon coupling. 

In Fig. H we have plotted the results for the dynamical 
structure factor, 



1 Im [D x (u 

7T 



- iS) + D 2 {uj + iS)] 



1 - 



exp(— /3ui) 



(13) 



where we have used the experimental resolution function 
[THz] 6 w 0.023 + 0.028o;/(27r) jl(J. The intensity of the 
experimental spectra |5j, also shown in Fig. [|, have been 
scaled; the (constant) background has been adjusted [ jlT| . 
The overall agreement between experiment and theory is 
satisfactory, although the hardening of the lower phonon 
mode is somewhat more pronounced in the experiment 
(6% vs. 1%). No experimental data for the upper mode 
were available for T = 16 K. In the inset a blowup of 
the central peak is given. It should be possible to re- 
solve the predicted central peak, which has a width of 
« 0.05 THz = 0.2 meV, by neutron scattering, testing 
thereby the theory. 
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FIG. 3. Theoretical (thick solid lines) and experimental 
(solid circles, |^|) results for the dynamical structure factor. 
The T 2 phonons are Peierls active and included in the theory, 
(XOX)i and (XOX)2 are other nearby phonons. The data for 
T = 296 K has been shifted. Inset: Blowup of the central 
peak at T = T S p = 14.1 K and T = 16 K (not broadened). 

The theory presented here is based on the RPA ap- 
proximation and we have shown it to yield reasonably 
good agreement with experiment, explaining the absence 
of a Peierls-active soft-phonon mode in CuGeC>3. From a 
theoretical point of view, one might question the applica- 
bility of RPA in the central-peak regime. A definite theo- 
retical resolution to this problem is not known at present, 
but one may note that the standard phenomenological 
theory for the central peak occurring in structural phase 
transitions has RPA- form Q. For the case of CuGeC>3 
we will test the applicability of RPA by comparing the 



prediction of RPA (with no fit parameter) for the pretran- 
sitional Peierls-fluctuations, given by the inverse lattice 
correlation length l/£, with the experimental results. 

The lattice correlation length is determined by the 
long-distance falloff, 



lim 

z — >oc 



dq 
2^ 



e iqz ReD q (0) 



z/c 



(14) 



where c = 2.94 A is the c-axis lattice constant of CuGc03 
and Dg{w) = Y^-y D qn (w). We have calculated l/£ from 
Eq. (g), using v s = (tt/2) J(l - 1.12a) @ (which en- 
ters Eq. (Q)), J = 156 K for the exchange integral and 
a = 0.24 for the frustration parameter [[l3| - |l5| . The re- 
sults for l/£ are presented in Fig. Q together with results 
for CuGeC>3 obtained by diffusive X-ray scattering |fl6| , 
which are consistent with neutron-scattering data and 
the absence of a soft phonon M. 
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FIG. 4. RPA and experimental results j^]] for the inverse 
lattice correlation length, l/£. The theory does not contain 
any free parameter. Inset: T 5 , ' AF ' (it, lj), as a function of 
lj/(2wT), as predicted by bosonization (solid line, Eq. (0)), 
and the neutron-scattering results [H], for T = 14.5 K (filled 
squares), T = 20 K (crosses) and T — 50 K (diamonds). 

Both experiment and theory show mean-field be- 
haviour, l/£ ~ y/T - T SP . The RPA-result agrees well 
with experiments. Above T — 19 K the lattice fluctua- 
tions have one-dimensional character Jl6| and the resid- 
ual difference between theory and experiment may be due 
to corrections to RPA. 

It is interesting to note, that Eq. (Q) for i s 
(for spin-rotational invariant Heisenberg chains) identi- 
cal with the bosonization result for the magnetic dynam- 
ical structure factor, (q, u>) (Tt]]. Quantum critical- 
ity implies T (n, iS) to be an universal function of 
w/(2wT), at least for small w. In the inset of Fig. ^ we 
have plotted the bosonization result for T M" 4 ^ (n, lo), 
together with the rescaled neutron-scattering results m 
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We observe that the experimental data approximately 
obey the scaling, though there is substantial scattering of 
the data for small ui/(2ttT), possibly influenced by non- 
critical contributions from the Peierls-fluctuations or by 
a crossover of the character of the magnetic excitations 



from ID to 2D near the Peierls transition 14|. It is also 
interesting to note, that the prediction for T S^ af \tt, lo) 
is independent of J and that the data for other ID 
Heisenberg antiferromagnets with very different values 
of the coupling J, like KCUF3 |18|, should fall onto the 



same universal curve presented in the inset of Fig. |4j. 
An experimental verification of quantum criticality for 
would imply also scale-invariance for Xn(w), 
since both coincide within bosonization ||[l7l. 

Generalization Until now we assumed spin-rotational 
invarianiance. Next we will show, that a central-peak 
regime occurs also in spin-Peierls transitions lacking spin- 
rotational invariance. An example is the spin-Peierls 
transition in a system of phonons coupled to an array of 
chains with Ising spins. This model was solved exactly 
by Pytte [fl9|| . It contains a parameter regime, where the 
transition is displacive and phonons do not become soft. 
In the oppositve limit, when the spin-chains are xy-like, 
the transition corresponds via the Jordan- Wigner trans- 
formation to the standard Peierls transition ^0|. Again 
one can show that soft-phonons occur in RPA, e.g. 
for T S p = J/10, only for Q a < 0.8 J. For Q a > 0.8 J the 
Peierls-active phonon does not become soft and a central 
peak arises at T$p- 

Discussion In this letter we have shown, that the 
RPA approach to the spin-Peierls transition includes 
both a soft-phonon and a central-peak regime. This re- 
sult is at first sight counterintuitive, as continuous lat- 
tice distortions below T$p are generally associated with 
a softening of the lattice above T$p- 

The eigenstates of the spin-phonon system evolve adi- 
abatically as a function of the spin-phonon coupling 
strength in the soft-phonon regime. In the central-peak 
regime a new magnetophonon appears at low frequencies 
and condenses at T$p, leading to the structural transi- 
tion and the formation of spin-singlets. This new collec- 
tive excitation is a superposition of a phonon with two 
magnons in a singlet state. The magnetophonon cou- 
ples to the phonon propagator and therefore shows up 
as a low-energy resonance in D q (uj), the central peak. 
The other resonance in D q (u>), at has the limit 
rim g5r _»o = Slo (tt) • Therefore, one usually regards w T 
to be the "true" phonon frequency. In terms of the eigen- 
states of the coupled spin-phonon system such a distinc- 
tion does not make sense. In the central-peak regime 
the spectral weight of D q (uj) is divided in between the 
"phonon-resonance" at u>„ and the soft magnetophonon. 

Conclusions The absence of a soft Peierls-active 
phonon mode in CuGeC>3 has been considered as a chal- 
lenge to theory. It has been argued [^2| that the Cross 
and Fischer theory is essentially incomplete, i.e. not ap- 



plicable to CuGe03. Here we point out, that the absence 
of soft phonons does actually find a natural explanation 
within the CF-approach. The calculated temperature- 
dependence of the phonon modes and that of the pre- 
transitional Peierls-fluctuations are in excellent agree- 
ment with experiment. A central peak of width 0.2 meV 
is predicted to appear at T$p- Finally we have pointed 
out, that a key ingredient of the theory, the quantum 
criticallity of Xn(w), can be tested, albeit indirectly, with 
neutron scattering through a test of the scale invariance 
of the magnetic dynamical structure factor, S^ AF ^ (tt , u>) . 
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